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After a short elementary introduction to the exact renormalization group for the effective 
action, I discuss a particular truncation of the hierarchy of flow equations that allows for the 
determination of the full momentum of the n-point functions. Applications are then briefly 
presented, to critical 0{N) models, to Bose-Einstein condensation and to finite temperature 
field theory. 



I. INTRODUCTION 

The exact, or non-perturbative, renormalization group (RG) [llsl] stands out as a very 
promising formalism to address non-perturbative problems, i.e., problems whose solution 
does not appear to be expressible as an expansion in some small parameter. It leads to 
exact flow equations which cannot be solved in general, but which offer the possibility for 
new approximation schemes. It has been applied successfully to a variety of problems, in 
condensed matter, particle or nuclear physics (for reviews with various points of view on the 
subject, see e.g. |6h11||). 

When only correlation functions at small momenta are needed, as is the case for instance 
in the calculation of critical exponents, a general approximation method to solve the infinite 
hierarchy of the flow equations has been developed Baa. This method is based on a 
derivative expansion of the effective action. However, in many situations, this is not enough: 
in order to calculate the quantities of physical interest, the knowledge of the full momentum 
dependence of the correlation functions is mandatory. 

The present paper deals with this issue, and summarizes work that has been done following 
the original suggestion by Blaizot, Mendez-Galain and Wschebor (BMW) [l^ . 1^ to obtain 
the momentum dependence of n-point functions from the flow equations. The strategy put 
forward in [2] is based on the fact that the RG flow at scale k involves the integration of 
fluctuations with momenta q < n, as insured by the presence of a cutoff function Rf^(q). 
Since this cutoff function also guarantees that the vertex functions are smooth functions 
of the momenta, these can be expanded in powers of jr? . The "leading order" (LO) of 
the approximation scheme consists in a truncation at the level of the flow equation for the 
two-point function ri^\ setting g' = in the vertices rjf^ and F^f^ that appear in this flow 
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equation. Doing so, and working in a constant external field, one can then express and 
tI^^ as derivatives of r|f ^ with respect to the background field, thereby closing the hierarchy 
of flow equations. Next-to- leading orders are defined by similarly truncating the hierarchy at 
the level of higher point functions. 

The price to pay is that the flow equations become differential equations with respect 
to a uniform background field, with integral kernels that involve the solution itself. These 
non-linear integro-differential equations are a priori difficult to solve. It is possible to do so, 
however, with a numerical effort comparable to that involved in solving the flow equations 
that result from the derivative expansion. We shall provide here examples of results obtained 
with this method in various contexts. 

The outline of the paper is as follows. In the Sect. [TI]we briefly recall some basic features 
of the exact renormalization group, and derive the flow equation for the effective action. We 
also digress on the use of the variational principle to obtain flow equations in the context of 
the hamiltonian formalism often used in dealing with non relativistic systems. Then, in Sect. 
mH we discuss approximation schemes, such as the local potential approximation, and the 
BMW approximation scheme. Connections with the formalism of the two-particle irreducible 
effective action are also discussed. In Sect. \IV\ we present three applications within scalar 
field theory with 0{N) symmetry in which the full momentum dependence of the 2-point 
function plays an essential role. The paper ends with a short conclusion. 

II. EXACT RENORMALIZATION GROUP FLOW EQUATIONS 

In this section we give an elementary introduction to the exact renormalization group, 
focussing on the effective action. For most of the discussion (except for the digression in 
subsection III Bp we shall restrict ourselves to the case of a scalar field theory with the classical 
action 



whose parameters m and u are defined at some "microscopic scale" A, to be specified. 

A. Flow equations for the n-point functions and their generating functional 

The strategy of the version of the renormalization group that we consider here is to 
build a family of theories indexed by a continuous parameter k, with the dimension of a 
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momentum, and such that fluctuations are smoothly taken into account as k is lowered from 
the microscopic scale A down to 0. In practice, this is achieved by adding to the original 
Euclidean action S a (non-local) term, quadratic in the fields, of the form 

The regulator, or cut-off function, R^{q^) is chosen so that: i) RK{q ^ k) ~ which 
effectively suppresses the contribution of the modes ip{q < k) by giving them a mass ~ k; ii) 
it vanishes rapidly for q ^ k, leaving the modes (p{q > k) unaffected. The explicit form of 
the cut-off function i?K(g) used in actual calculations will be specified later. 

We consider then the "deformed" (non local) field theory with action 3^ = 8 + AS^, and 
the functional integral that yields the generating functional of Green's functions 

ZM = I V^e-^'^+f^^, (3) 

with J j(p = J dxj{x)ip[x). The expectation of the field in the presence of the source j{x) is 
given by 

—-f = {^{x)=cl>{x). 4 
6j[x) 

Similarly 

d^lnZ, = -^l dMQ){^iQM-Q)), (5) 



where we used here the fact that all the dependence on k is contained in the regulator term 
AS'k[93]. Note that the expectation values in Eqs. dH and ([5]) are taken in the presence of 
the regulator, and depend therefore on k. At this point, we define the effective action for the 
deformed theory, T^ [</>], through a Legendre transform 

r«M+lnZ.[j]= ^^=J- (6) 



The functional [(j)] is the generating functional of the one- line irreducible n-point functions 
of the deformed theory. Taking into account that, at fixed (p, j depends on k, one easily 
obtains the fiow of ^^.[(j)] 

dTni'P] = 11 dMq){¥^iqM-q))- (7) 
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It is in fact convenient to redefine T^ by subtracting from it A^^. This subtraction has the 
advantage to make coincide with the classical action S at the microscopic scale A (rather 
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than to make it coincide with Sa = S + AS a, as would be the case for the definition (^). 
Furthermore, it modifies the flow equation in such a way that only the field fluctuations are 
involved, 
function 



hat is, once the subtraction is made, the flow is driven by the connected 2-point 

fl: 



d^^M = \j dM^) G.{q, -q; (8) 

where 

G{q, -q; 0) = {^{qM-q))c = {^{qM-q)) " Hq)H-Q), (9) 
is the full propagator in the presence of the background field 0. Using well-known properties of 

(2) 

the Legendre transform, one can relate G{q, —q; (jj) to [(/>], the second functional derivative 
of Tj^cf)] w.r.t. 4>: 

G-iM=rl2)[0]+i?,. (10) 

The initial conditions on the flow equation ([8]) are specified at the microscopic scale k = A 
mentioned above. This scale A is the scale where fluctuations are damped by A5'k, so that 
rf5=A[</'] ~ S[(j)\. It should be emphasized however that this issue of the initial condition at 
the scale A is a subtle one. It is intimately related to that of renormalization and ultraviolet 
divergences, issues that will not be discussed here. The effective action of the original scalar 
field theory is obtained as the solution of Eq. ([8]) for k — )• 0, at which point R^^^q^) vanishes. 
Whether at that point the result obtained for T[(j)\ is independent on the choice of R,^, that is 
on the path joining the classical action to the full effective action, is another subtle question, 
whose answer depends on the approximation made. We shall comment on it further when we 
discuss specific approximations. 

By taking successive functional derivatives of eq. ([8]) with respect to (p, and then letting 
the field be constant, one gets the flow equation for the n-point functions 



^ (xi , • • • , Xn , 



(11) 

(x)=ct> 



6(t){xi) . . . 6(j){Xr, 

in a constant background field. Since the background is constant, these functions are invariant 
under translations of the coordinates, and it is convenient to factor out of the definition of 
their Fourier transform the (^-function that expresses the conservation of the total momentum. 
Thus, with an obvious abuse of notation, we define the n-point functions T^^\pi, . . . ,pn;ct>) 
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as: 



= j (fxi... j <fxne'^UP^^^T^^){xi,--- ,Xn-A)- (12) 

We use here the convention of incoming momenta, and it is understood that in 
ri"'* (pi , . . . , p„ ; (;/)) the sum of ah momenta vanishes, so that T^k\pi, . . . ,pn',(l)) is actuahy 
a function of n — 1 momentum variables. We shall often use the simplified notation T\_ (p; (f)) 
for the function (p, — p; (j)). 

When (j) is constant, the functional itself reduces, to within a volume factor to the 
effective potential 14(0): 

[(/)] = riVK [</>], (/) constant. (13) 

The flow equation for the effective potential V^, follows from that of the effective action F^, 
Eq. dH]), when restricted to constant (p. It reads 

i^dMi') = \j ^ Gn{q, <A), (14) 

where 

G-\q,<P)^r(^\q,cl>)+R,iq). (15) 

By taking two derivatives of eq. ([8]) with respect to <f), and then letting the field be constant, 
one obtains the equation for the 2-point function: 

9«F(2)(p,</>) = I ^dMq)GUq,cp) 

rL') {p, q, -p - q; cl))G^{q + p, ^)r'i\-p, p + q, -q; 4>) - ^FI'^) {p, -p, q, -q; 0) 



(16) 



A diagrammatic illustration of this equation is given in Fig. [TJ The flow equations for the 
n-point functions, such as Eq. (I16p . constitute an infinite tower of coupled equations. The 
coupling between equations occurs in two ways, upwards and downwards. Upwards: typically 
the equation for F*-"-* involves F^""*"^^ and F^""^^-*. Downwards: all the flow equations involve 
F*^^) which is coupled successively to all the equations above it. There exist also consistency 
conditions that we illustrate here with F^^^ and the effective potential: T^'^\q = 0) = d'^V/dcfP' . 
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FIG. 1: The two diagrams contributing to the flow of the 2-point function, Eq. p^ . The internal lines 
represent dressed propagators, Gk- The cross represents an insertion of d^Rk- The vertices denoted 
by black dots are rlf^ (left) and V^^^ (right). 

Thus, r(2)(g = 0) can be calculated by taking the second derivative with respect to (f) of the 
flow equation for the effective potential: 



X < 



{2-kY 



(17) 

Alternatively, it can be obtained from Eq. ()16p where we set p = 0. The two calculations 
agree owing to the fact that 

rf((?;(A) .4^, ,,,, 92rl=^^((Z;<A) 



r'i\q,-q,0;< 



r<^^\q,-q,0,0;. 



(18) 



We shall consider later approximation schemes that treat at different levels of accuracy the 
equations for various n-point functions. This may lead to violations of such consistency 
conditions if one is not careful. 

Relations such as (jl8p between n-point functions with some vanishing external momenta 
and derivatives of n-point functions of lower rank with respect to a constant background field 
play an important role in the approximation scheme presented below. In this context, we 
simply note here that the effective potential may be viewed as the generating functional of 
the n-point functions with vanishing external momenta. That is (for constant (p) 



rl")(o,--- ,0) 



(19) 



(2) 

In the same way, we may regard [p] (j)) as the generating functional of n-point functions 



with one non-vanishing momentum 



rW(p,-p,o,--- ,0) 



d^-^T^n\p;ct>) 



kn-2 



(20) 
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B. Flow equation from the variational principle 

The previous derivation of the flow equation did not refer to the interpretation of the 
various manipulations that were done in terms of "coarse graining" or "ehmination of degrees 
of freedom" , that are essential aspects of the renormalization group. These are actually hidden 
in the specific functional form of the regulator Rf^(q), but nothing in the derivation of the 
flow equation actually depends on the choice of the regulator (asides from being quadratic in 
the fields, but even that can be generalized). Thus, one may view the fiow equation as merely 
a tool to continuously go from a "simple" action (the classical action) to the full effective 
action, by tuning an appropriate contribution to the action. This point of view opens new 
perspectives and suggests that the strategy of fiow equations may have a larger flexibility 
than is presently exploited. As an illustration, we shall here indicate how the variational 
principle can be used to establish flow equations in terms of hamiltonians and wave-functions 
for non-relativistic systems. 

The strategy we shall follow is identical to that of the previous section. The system is 
described by a hamiltonian H, to which is added an "external held" Ai^^ that plays the 
role of the regulator. We assume that AHf^^o = 0, and that H\ = H + IS.H\ is somehow 
"simple", so that for instance one can determine its ground state. Then the flow equation 
will take us from the ground state of this simple hamiltonian H\ to the ground state of the 
hamiltonian H. 

In this context, the initial condition of the flow is not restricted, as it is in most field 
theoretical applications to be some "classical action" . Thus, it may be advantageous to depart 
from the standard approaches that would take Hj^ to be the hamiltonian of independent 
particles, or quasi-particles, and include part of the effects of the interactions into H\. In 
fact such a strategy was recently implemented by Machado and Dupuis, in their lattice 
renomalization group [l^, in which they separate the on-site physics form the long range 
correlations, only the latter being treated by the renormalization group. Much earlier, Parola 
and Reatto have developed a theory of liquids in which the short range part of the interaction 
is treated (almost) exactly with hard sphere systems, while the long range correlations are 
treated by a fiow equation [l^, 3] that actually resembles closely Eq. ([8]). 

The fiow equation for the effective action, Eq. ([8]), is general and could be used in the 
non-relativistic many-body context. However, to bring a new perspective to the discussion, 
we shall go through an extremely simple derivation of a fiow equation that relies on the 



variational principle [ITy. For simplicity, we shall assume that Aff^ is a one body operator, 
i.e., a quadratic form of creation and destruction operators, that depends on the continuous 
parameter k: 

AH, = Y,RKiQ)aia^, (21) 
q 

with Oq and Oq creation and destruction operators (of fermions or bosons). We wish to 
calculate the ground state energy Ef^ of the hamiltonian H,^ = H + AH^. We call \^k) the 
corresponding eigenstate. In fact we do not need the exact eigenstate, it is enough that 
be determined from the variational principle 

6{^^\H + AH^\^^)\ji^=0, (22) 

where it is understood that the regulator i?K remains fixed in the variation. The flow equation 
for E'k follows then immediately 

dJ^ = Y,dMq){'^.\a\a^\^^), (23) 
q 

where (^«;|fflqOq|^K) can be seen as an occupation factor. The equation above is the analog 
of Eq. ([5]). It just describes the change in the ground state energy under the change of the 
parameters of the hamiltonian. Of course, when some creation or destruction operators ac- 
quire expectation values (as in the case of Bose-Einstein condensation), it may be convenient 
to work with the connected part of (^'fclapapl^K;), as we did earlier. Note that the flow of 
the energy = (^'kI^^I^'k) reads 

d,E^ = -Y,RMdn{^n\alag\^,). (24) 
q 

This approach is useful if we can determine (approximately) the wave function, or at least 
the occupation factors, as a function of k. It was applied in the study of the BCS-BEC 
crossover [17]: in this case the BCS wave function leads to a reasonably accurate estimate of 
the occupation factors and their dependence on k, and this extremely simple scheme allowed 
us to reproduce qualitatively (and even semi quantitatively) the results of more complete 
renormalization group studies, such as those carried out in Ref. [3] • 



III. APPROXIMATION SCHEMES 



The flow equations for the n-point functions that have been presented in the previous 
section are exact. Their solution requires, in general, approximations. It is precisely one of 
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the virtues of the formulation of field theory based on the exact renormalization group to 
suggest approximations that are not easily derived in other, more conventional approaches. 
Particularly interesting are the approximation schemes for the effective action itself, that is, 
approximations that apply to the whole set of n-point functions at once. The approximation 
schemes to the discussed in this section have this property. 



A. The local potential approximation 

The local potential approximation (LPA) is the simplest of such approximations. It con- 
sists in assuming that for all values of k the effective action takes the form 

TM = I d'x |i (90)2 + K(</')} . (25) 



Within this approximation, the 2-point function, obtained by differentiating Tfi[(p] in Eq. (|25p 
twice with respect to (p, and letting (p be constant, is of the form (p = (p'^/2) 

r'iHq; p) = + mlip) ml{p)^^. (26) 

The corresponding propagator is simply a massive propagator, with a (}!)-dependent mass 
niKip)- "T^G flow equation for the effective potential becomes a closed equation that we 
write as follows 

K9,K(p) = ^/i(At,p), (27) 

where Ii{k, p) is given by 

In{K,p) = Jn{p = 0,K,p), J^{p, p) = J kOMq) G .{p + q)G':-\q) . (28) 

Note that the form (j27p of the equation for the effective potential is general: it would yield 
the exact effective potential if the propagator used to calculate Ii was the exact propagator, 
instead of the LPA propagator. Note also that, within the LPA, the equation for the 2-point 
function at vanishing external momentum, given by Eq. (jlZp . reduces to a closed equation 

Kd^Tl,>{0;p) = l3{K,p) \ — I --l2{K,p) , (29) 

(2) 

since the derivatives of Tk {q] p) with respect to (j) are identical to derivatives of the effective 
potential, or, equivalently, to derivatives of Fj^ ^(0; p). 
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The LPA has 
is quite good 



aeen widely used, and given its simpUcity, the quahty of the results obtained 
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23l |. Note that the LPA, by construction, yields n-point functions that 
have no intrinsic dependence on the external momenta: the generating functional of these 
n-point functions is given by Eq. (j25p taken at k = 0, so that, for n > 2, the LPA n- 
point functions are just derivatives of the effective potential (see Eq. (I19p ). The LPA can be 
improved through an expansion in gradients of the fields, usually referred to as the derivative 
expansion (DE) 6|, |8|, |20|. Within the derivative expansion, the lowest n-point functions 
have polynomial dependence on the momenta (the degree of the polynomial corresponding 
to the order of the expansion in the derivatives). It follows in particular that the derivative 
expansion does not directly describe the anomalous behavior of n-point functions at small 
momenta (anomalous dimensions can be recovered for instance from the scale dependence of 
the field normalization). In order to capture the full momentum dependence of the n-point 
function, a better truncation scheme is necessary: the BMW scheme [12I ]. to be discussed 



next, achieves this goal. 



B. The approximation BMW-LO 

The BMW approximation scheme relies on two observations. First, the presence of the cut- 
off function Rk{q) insures that the n-point functions T^f^\pi) remain regular functions of the 
external momenta pi as pi — t- 0; besides, it limits the internal momentum q in equations such 
as Eq. (|16p to q ^ k. In line with this observation, the approximation consists in neglecting 
the (/-dependence of the vertex functions in the r.h.s. of the flow equations (e.g. set q = in 
r^^) and r(^) in Eq. ()16p ). while keeping the full dependence on the external momenta pi. The 
second observation is that, for uniform fields, ^^i^^^\pi, ■ ■ ■ ,Pm, 0, (/)) = c^^r^™"^ (pi, • • • ,Pm, </*) 
(a relation that we have already used earlier, see Eq. (llSp ). This enables one to close the 
hierarchy of equations at some finite order. The order m of the scheme consists in keeping 
the full momentum dependence of all the n-point functions up to ri™^ , and expressing p^,'""''^) 
and p(,'"+^) as derivatives of rfj"^ with respect to (p, after setting to zero the loop momenta 
that fiows through rl™"^^^ and p(,™+^) [n the equation for ri^\ The accuracy of the scheme 
depends of course of the rank m at which one operates the truncation, but obviously the 
implementation becomes increasingly complicated as m grows. It is therefore gratifying that 
accurate results can be obtained with lowest order truncations. 

The approximations obtained by truncating the hierarchy at the lowest level, i.e., m = 
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is identical to the local potential approximation discussed in the previous subsection. The 
next order of the approximation scheme, which we shall in fact refer to as the leading order 
(LO) of the BMW method, consists in a truncation at the level m = 2, i.e., at the level of 
the flow equation of the 2-point function. The resulting equation is the closed equation 

where J3{p,K,p) and I2{k, p) are obtained from the general definitions in Eq. ()28p . This 
equation may be viewed as a generalization of Eq. (I29p that takes into account the full 
momentum dependence of the 2-point function. What the BMW approximation achieves is 
the possibility to factorize the vertices and take them out of the integrals. The momentum 
dependence that remains within the 3 and 4-point vertices is that of the 2-point function 
itself (see Eq. (j20p ). Interestingly, a very similar set of equations have been obtained much 
earlier in the context of the theory of liquids ^] , but these were unknown to the authors of 
Ref. 0. 

At this point we recall the consistency condition discussed in connection with Eq. (|17p : 
the flow equation for {p = 0,p) can be obtained by taking the second derivative of the 
flow equation for the effective potential with respect to the background field. The resulting 
equation, Eq. (fTTl) . does not coincide here with Eq. (|30]l in which we set p = 0. This is 
because Eq. (j30p results from the BMW truncation, which is not implemented in Eq. ()17p . 
and needs not be. This results in a mismatch between the two ways of calculating T^'^\p = 
0), whose origin can be traced back to the fact that the 2-point function and the effective 
potential are not determined with the same accuracy: loosely speaking, the 2-point function 
calculated from Eq. ()30p is accurate to order one-loop, while the effective potential obtained 
from Eq. (I27p . with the BMW propagator, is accurate to order two-loop. In order to properly 
deal with this feature, we treat separately the zero momentum (p = 0) and the non-zero 
momentum (p ^ 0) sectors, and write 

r^^\p,p)^p^ + A^{p,p)+ml{p), (31) 

where A^{p = 0,p) = 0. Now, mf.{p) is obtained by solving the equation for the effective 
potential, while the equation for Ak(p, />) can be easily deduced from that for , i.e., from 
Eq. (i30]) by subtracting the corresponding equation that holds for p = 0. It reads 

dtA^{p,p) = 2pJ3{p,K,p) [u^,{p) + A'^{p, p)]"^ - 2pl3{K, p) ul{p) 

- ^l2(K,p) [A;(p,p) + 2pA'^(p,p)] , (32) 
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where the symbol ' denotes the derivative with respect to p, and Ufj_(p) = dmf.{p)/dp. 

This equation (I32p . together with that for the effective potential, Eq. (I27p . and that for 
the propagator 

G-\q,p) = q^ + A^{q,p)+ml{p)+R^{q), (33) 

constitute a closed system of equations for [p, p) which can be solved with the initial 
condition (p, p) = + m? + up. 



C. Relation to the 2PI formalism 



Because the exact RG formalism that we are using puts emphasis on the propa gato r, it 



is natural to ask for the connection with the 2PI (2-particle irreducible) formalism [244261] . 
Let us recall that the central quantity in this formalism is *1>[G], the sum of the two-particle- 
irreducible "skeleton" diagrams, a functional of the full propagator G. From ^[G] one obtains 
the self-energy by functional differentiation (to within factors (27r)^): 

This relation, together with Dyson's equation: 

G-\p)=p^ + + T.{p), (35) 

defines the physical propagator and self-energy in a self-consistent way. We shall refer to 
Eq. (p5]) . with S[G] given by Eq. as the "gap equation". A further differentiation of 

^[G] with respect to G yields the two-particle- irreducible kernel 

of a Bethe-Salpeter type equation 

F(^)(5,p) = X{q,p) - jT^^\q,l)G\l)X{l,p), (37) 

that allows the calculation of the four-point function T^^\q,p) = T^^\q, —q,p, —p)'- the quan- 
tity I{q,p) is the two-particle-irreducible contribution to T^'^\q,p) in one particular channel. 
If all skeletons are kept in these relations are exact. A ^-derivable approximation [25] is 
obtained by selecting a class of skeletons in $ and calculating S and F^^^ from the equations 
above. 
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The 2PI formalism provides a set of functional relations among the n-point functions that 

n 

can be used to define a truncation of the flow equations [27]. Consider indeed the equation 
for the 2-point function for a vanishing background field. It reads: 

9«r(2)(p)= - ^l^dMq)GUq)r(tHq,p). (38) 

A possible truncation consists in using for rif^(g,p) in the right-hand-side of this equation, 
the extension of the relation (j37p to the deformed theory, namely 



r(f)(g,p)=X«(g,p) - ^J^Ti^\q,l)Glil)Ul,p), (39) 

where the subscript k on means that the functional derivative defining the kernel I (see 
Eq. dMD) is to be evaluated for G = G^, with G'^ = T^^^ + R ft. Since X,j is a functional of 
the 2-point function, the system of equations (j38ti39p is indeed closed. 

One nice feature of this truncation scheme is that it is systematically improvable, by 
adding more skeletons to ^: if all skeletons are included, the solution of the coupled system 



of equations (|38ll39p provides the exact 2-point function as well as the exact 4-point function 
for a particular configuration of the external momenta. A second attractive feature is that it 
preserves the property of the flow of being a total derivative with respect to the parameter 
K. It is indeed not difficult to show, using Eqs. (f39l) and (f38l) . that 



a,r(2)(p) = 5,s«(p), (40) 

where = with given by Eq. This is a unique property of this truncation, 

that is not shared by most other popular truncations of the exact RG (with the noticeable 

exceptions of the perturbative expansion, and the large approximation, see e.g. [ssl- A 

similar property of the flow equation was also obtained in the off-equilibrium context in 
I I 

Ref. [30]. The BMW truncation does not respect this property). Since it is a total derivative, 
the flow can be easily integrated to yield the gap equation whose solution is equivalent to a 
resummation of all the Feynman diagrams that are generated from the skeletons that are kept 
in the approximation considered. And because the solution of the gap equation corresponds 
to an exact resummation of selected Feynman diagrams, at the end of the flow where k = 
and the regulator vanishes, the final result is rigorously independent of the choice of the 
regulator. 

All these properties of the 2PI truncation may look at first disappointing from the point 
of view of the flow equations: indeed, all what the flow does in this particular truncation is 
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solving the 2PI equations! However, there is certainly interest in establishing direct connec- 
tions between non trivial non-perturbative approximations. In particular, because the 2PI 
truncations lead to flow equations that are exact derivatives, they could be used to test other 
approximations. Besides, from the point of view of the 2PI formalism, there is a practical 
advantage in reformulating the gap equation as a flow equation: this is because initial value 
problems are in general easier to solve than non linear gap equations. Furthermore, regarding 
the 2PI equations as flow equations shed a new light on the renormalization of <I>-derivable 



approximations 



2S 



IV. APPLICATIONS 

We turn now to a few specific applications of the non perturbative renormalization group 
using the BMW truncation scheme. Most of these applications concern scalar field theories 
with an 0{N) symmetric action of the generic form 

S = J d'x}^-{dMx)f + ^^Hx) + ^{y.\x)fj, ^' = Y^^^a^,. (41) 

The first application covered in this section may be seen as a "classic" application within 
the context of the renormalization group: it concerns the critical 0{N) models, and the 
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32| . The second application 



calculation of critical exponents and the scaling functions 
refers to the shift of the transition temperature of the Bose-Einstein condensation in a dilute 
Bose gas: this again involve the 0{N) field theory (with N = 2), but the relevant quantity to 
be calculated is sensitive to the whole momentum range of the 2-point function, and not only 
to the critical momentum region 



33l ]. The last application concerns the thermodynamics 
of quantum fields. It shows how the exact renormalization group allows one to circumvent 
the specific difficulties of perturbation theory in such systems. The physics motivation is 
the physics of the quark-gluon plasma, but some of the essential difficulties of perturbation 
theory in QCD at finite temperature are well illustrated again by the scalar field 



34 



35| 



Before turning to these applications, we need to comment on the dependence of the reg- 
ulator in the practical calculations. We shall use in our calculations a cut-off function of the 
generic form 

R^iq) = Z^K^iq), q^^, (42) 

K 

where is a function of k only. This factor reflects the finite change in normalization of 
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the field between the scale A and the scale k. It can be defined as 



(43) 
po;po 



with pq and po a priori arbitrary. In practice, one usually chooses po = and po the value of 
p at the minimum of the effective potential. The factor enters the scaling to dimensionless 
variables used for the numerical solution in the critical region. 

The results of the calculations presented below were obtained with an exponential regu- 
lator 

r{q) = (44) 
e'J — 1 

where a is a free parameter. As we shall see, physical quantities exhibit a small dependence 
on a. Since in the absence of any approximation, they would be strictly independent of the 
cut-off function, a study of this spurious dependence provides an indication of the quality of 
the approximation P, [fi^] • 



A. Critical 0{N) models 



The equations of the BMW method have been solved first with additional approximations 



in Refs. 



19[ , but the results presented here were obtained by solving numerically the nonlinear 



integro-partial-differential equations (j32p and (j27p without any further approximation 3l|. 



32] • Here, we report only some 



The numerical techniques used are described in Refs. [31| and 
selected results, and comment on important aspects of the approximation that are needed in 
order to gauge the quality of these results. 

Table U contains results for the critical exponents rj and u, in dimension d = 3 and for 
various values of N, together with some of the best estimates available in the literature. 
To appreciate the quality of these results, let us recall that these depend a priori on the 
parameters a, po, and pQ (see Eqs. ([l2]) and (H3]l ). (They also depend on the functional 
form of the cut-off function, but no systematic exploration has been performed to study this 
dependence.) In all cases studied, we find the dependence on pQ and po to be much smaller 
than that on a, so that only the latter needs to be considered. As a function of a, physical 
quantities typically exhibit a single extremum a* located near a = 2. Moreover, we find this 
extremum generally pointing towards the best numerical estimates. Since, in the absence of 
any approximation there should be no dependence on a, we regard these extremum values. 
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TABLE I: Critical exponents for the 0{N) models for d = 3. 



N 


BMW 


Resummed pert. exp. 




MC and high-temp. 


series 




V 


V 


r\ 


V 


Ref. 


r/ 


V 


Ref. 





0.034 


0.589 


0.0284(25) 


0.5882(11) 


m 


0.030(3) 


0.5872(5) 


[38] 


1 


0.039 


0.632 


0.0335(25) 


0.6304(13) 


m 


0.0368(2) 


0.6302(1) 


[37] 


2 


0.041 


0.674 


0.0354(25) 


0.6703(15) 


[36] 


0.0381(2) 


0.6717(1) 


m 


3 


0.040 


0.715 


0.0355(25) 


0.7073(35) 


m 


0.0375(5) 


0.7112(5) 


mm 


4 


0.038 


0.754 


0.035(4) 


0.741(6) 


[36] 


0.0365(10) 


0.749(2) 




10 


0.022 


0.889 


0.024 


0.859 


[42] 









being locally independent of a, as our best values, adopting in doing so a strategy often 
referred to as the principle of minimal sensitivity (PMS) 4J]. Our numbers are then all 
given for the PMS values d* of the regulator parameter, and the digits quoted are those 
which remain stable when a varies in the range [a* — \-,ol* + ^]. The quality of the values 
obtained for the critical exponents is obvious: For all N the results for v agree with the 
best estimates to within less than a percent; as for the values of r\ it is typically at the same 
distance from the Monte-Carlo and temperature series estimates as resummed perturbative 
data. For N = 100, we find r] = 0.0023, and v = 0.990, which compare well to the values 
r] = 0.0027 and v = 0.989 obtained in the 1/A^ expansion 4511 . Our numbers also compare 
favorably with those obtained at order 9^ in the DE scheme 2C ] . 

The two-dimensional case, for which exact results exist, provides an even more stringent 
test of the BMW scheme. We have results at the moment only for the Ising model, i.e., for 

= 1, which exhibits a standard critical behavior va. d = 2. The perturbative method that 
works well in d = 3 fails here: for instance, the fixed-dimension expansion that provides the 
best results in d = 3, yields, in d = 2 and at five loops, r] = 0.145(14) 40], in contradiction 
with the exact value rj = We find instead rj = 0.254, v = 1.00. Note however that 
no systematic study of the dependence of this result on the regulator parameter has been 
performed yet. 

The BMW scheme yields the complete momentum dependence of the 2-point function. All 
expected features of F^, (p) at criticality are observed: In the infrared (IR) regime k <^ p <^ u, 
Tli\p,0) ~ p^~'^k'^, and this IR behavior of tI^\p) can be used to extract the value of 
rj] the value obtained directly from the momentum dependence of F^ (p) is in excellent 
agreement with that deduced from the K-dependence of the field normalization factor Z^. 
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FIG. 2: The scaling function g+{x) as a function of x = q^. The experimental data are from Ref. 

The ultraviolet (UV) regime n,u <^ p <^ A exists if u is sufficiently small; this regime can 
be studied perturbatively and one finds that, in leading order, t\^\p,0)) ~ u'^log{p/u). The 
present approximation reproduces this logarithmic behavior with, however, a prefactor 8% 
larger than the two-loop result. Note that the complete two-loop behavior can be recovered 



32| 



by a simple improvement of the BMW scheme 

The momentum dependence of the 2-point function can be further tested by analyzing 
the so-called scaling functions 



G±iq) = X9±{q0: 



(45) 



where = r(^^(g = 0), and ^ is the correlation length, which diverges close to criticality 
with the I' critical exponent. Here it refers to the two phases, above and below the critical 
point, respectively. The functions g± are universal. We consider here the case = 1. The 
scaling function = g^) has been calculated in the BMW approximation for different 

values of the correlation length. It is plotted in Fig. [21 where we can see that the scaling is 
perfectly reproduced. The agreement with the experimental data is also excellent, although 
this is somewhat delusive since, in this range of momenta, the physics is dominated by 
mean field effects, and the scaling function deviates only slightly from that predicted by the 
Ornstein-Zernicke approximation). A critical study is presented in Ref. 32l |. 
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B. Temperature of Bose-Einstein condensation 

A quantity particularly sensitive to the UV-IR crossover region is the shift, due to in- 



43]. The 



teractions, of the critical temperature of the dilute, weakly interacting, Bose gas 
hamiltonian describing such a system is typically of the form 

H = j d'v l^t(r) ^(r) + |v.t(r)Vt(r)V'(r)V'(r)| , (46) 

where g = ATra/m, with a the s-wave scattering length. (We are ignoring here a subtlety 
related to the ultraviolet divergences that are generated in calculations with a contact inter- 
action and that require the introduction of an ultraviolet cut-off. This can be implemented in 
a standard fashion, but it plays no role in the present discussion.) The effective hamiltonian 
(j46p provides an accurate description of dilute systems, when the scattering length is small 
compared to the interparticle distance, i.e., an^/^ <ti 1. In the vicinity of the condensation, 
where n^/^A ~ 1, with A = \J~^ the thermal wavelength, the diluteness condition reads 
a/A < 1. 

It has been shown that the shift ATc of the Bose-Einstein condensation temperature is 
linear in an^/^ 48 1: 



-'^c ^ c 

Here is the condensation temperature of the ideal gas, defined by the condition nA^ = 
C(3/2) ~ 2.612, and Tc the transition temperature of the interacting system at the same 
density. This result is non trivial: although the shift is proportional to a, and hence is small 
if a is small, the result (j47|) cannot be obtained form perturbation theory. It is obtained 
through the following chain of arguments. In the limit of small coupling the shift of the 
critical temperature may be obtained from the shift of the critical density, = — |^^- 
The latter quantity is easier to calculate, and is dominated by the contribution of the zero 
Matsubara frequency part, -001 of the bosonic field, whose dynamics is governed by a three- 
dimensional classical field theory with 0(2) symmetry. To make the connection with the 



notation used in the rest of this section, we set Vo = Vn^T{ipi + iip2), with ipi and ip2 two 
real fields. The action for these real fields is given by Eq. (I41jl . with the parameters r and u 
related to the parameters of the hamiltonian (j46p by 



-2mTa, tt = 967r^-^. (48) 
A^ 
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The resulting density shift is then given by the change in the fluctuation of this classical field, 
with the coefficient c in Eq. (I47p given by 

256.3 A(^;) 

in the limit ti — t- (and for N = 2). 

The calculation of A((/9?) is difficult for the reason already mentioned: although the cou- 
pling constant u can be arbitrarily small, perturbation theory cannot be used because of 
the infrared divergences of the critical three-dimensional field theory. The best numerical 
estimates for A{ip'^), and hence for c, are those which have been obtained using the lattice 



technique by two groups, with the results: c = 1.32 it 0.02 



5l| and c = 1.29 ± 0.05 [52]. The 



availability of these results has turned the calculation of c into a testing ground for other 
non perturbative methods: expansion in 1/iv 0,0, optimized perturbation theory 



55 



53|. 



resummed perturbative calculations to high loop orders 

To understand better the origin of the difficulty of the calculation of c, as well as the 
linearity in a of ATc, let us write A(93?) as the following integral 



AM 

N 



1 



1 

'2^ 



dp 
p 



P 



p. 



p^ + T,{p) 



(50) 



(27r)3 Vp2 + S(p) p2 

where S(p) is the self-energy at criticality, i.e., S(0) = 0. In eq. (j50p . the term within 
the square brackets is, to a very good approximation (when p/u is small enough), equal 
to Ti{p)/p, a function that is peaked in the region of intermediate momenta between the 
critical region and the high momentum perturbative region. Thus, the difficulty in getting a 
precise evaluation of the integral (j50|) is that it requires an accurate determination of S(p) 
in a large region of momenta including the crossover region between two different physical 



regimes 



4S 



57[. In that sense, the calculation of c can be viewed as a very stringent test of 



the approximation scheme, and in fact, it is in order to obtain a reliable estimate of c that 
the BMW scheme was initially developed. 

To see the origin of the linear relation between Aric and a, we note first that the action 
([T]) contains a single dimensionfull parameter, u, r being adjusted for any given u to be at 
criticality. In fact the effective three dimensional theory is ultraviolet divergent, so there is a 
priori another parameter, an ultraviolet cut-off A ~ 1/A, with A the thermal wavelength. It 
follows then from dimensional analysis that S(p) can be written as $](p = xu) = u'^a{x, u/A). 
Now, the diagrams involved in the calculation of S (at criticality) are ultraviolet convergent, 
so that the infinite cut-off limit can be taken. Note however that this requires that all 
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TABLE II: Coefficient c for the 0{N) models. 



N 


BMW 


Resummed pert 


th. 


Monte-Carlo 




c 


c 


Ref. 


c Ref. 


1 


1.15 


1.07(10) 


[53j 


1.09(9) [54] 


2 


1.37 


1.27(10) 


[53] 


1.32(2) [51] 


3 


1.50 


1.43(11) 


[53] 




4 


1.63 


1.54(11) 


[53] 


1.6(1) [54] 



momenta involved in the various integrations are small in comparison with A or, in other 
words, that the integrands are negligibly small for momenta k ~ A^^. Only then can we 
ignore the effects of non vanishing Matsubara frequencies, and finite cut-off effects. This 
implies that uX ~ a/X is sufficiently small. In this region of validity of the classical field 
approximation, that is, for small enough u, a{x,u/A) becomes a universal function a{x), 
independent of u, and 



Nu 



1 

2^ 



dx 



a{x) 



x"^ + (t{x) ' 



(51) 



showing that the change in the critical density is indeed linear in u, and hence in a. 

Table [TI] contains our results for c together with some of the best estimates available in 
the literature. As was the case for the critical exponents in Table [H our numbers are all given 
for the PMS values a* of the regulator parameter, and the digits quoted are those which 
remain stable when a varies in the range [a* — ^, q* + ^]. For all N values where six-loop 
resummed calculations exist, our results for c are within the error bars. For N = 100, we 
find c = 2.36, which compares well to the exact large N value c ~ 2.33 49|]. Our estimates 



for c are also coirroarable to those obtained from an approximation specifically designed for 



this quantity 



13|, 



60|. See also 58| for an earlier estimate of c using a different truncation of 



the RG equations, and also 



59|. 



C. Thermodynamics of quantum fields 



The last application that we shall consider concerns the thermodynamics of quantum 
fields at high temperature. This is motivated by the study of the quark-gluon plasma. In 
such a system one could expect a priori perturbation theory to yield accurate results because 
the asymptotic freedom of Quantum Chromodynamics makes the effective coupling small at 
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high temperature. However, strict perturbation theory does not work: it exhibits indeed 
very poor convergence properties, even in a range of values of the couphng constant where 
reasonable results are obtained at T = 0. This difference of behavior of perturbation theory 
at zero and finite temperature can be understood from the fact that, at finite temperature, 
the actual expansion parameter involves both the coupling constant and the magnitude of 
thermal fluctuations (for a recent review, see |61|]; see also [62])- In that respect, the problem 
is not specific to QCD: Similar poor convergence behavior appears also in the simpler scalar 
field theory 63|, and has also been observed in the case of large- ip^ theory 6j]. We focus 



351]. 



here on the case of scalar field with a g^^'^ interaction (i.e., cp' = n/24) 

Let us first recall how the effect of the interactions at a given scale depends on the 
magnitude of the relevant thermal fluctuations at that scale (and in some cases at a different 
scale as well) . The thermal fluctuations of the field are given by the following integral 

When we perform a perturbative calculation, we assume that the "kinetic energy" ~ {{d(p)'^) 
is large compared to the "potential energy" ~ g^{ip'^). Obviously, this comparison depends 
on the strength of the coupling, but also on the typical wavelength, or momentum, of the 
fluctuations. To make things more precise, let us observe that the integral (I52p is dominated 
by the largest values of k (in the absence of the statistical factor it would be quadratically 
divergent). One may then calculate the integral with an upper cut-off k and refer to the 
corresponding value as to "the contribution of the fluctuations at scale k", and denote it 
by {f'^)K- In the same spirit, we may approximate the kinetic energy of modes at scale k 
as ((5(/7)2)k Assuming furthermore that ~ one gets the expansion 

parameter (ratio of potential energy ~ g^{'f'^)'^ to kinetic energy ~ ^^{lp^),^) 

7. = (53) 

Let us then examine the values of this parameter for several characteristics momenta. The 
fluctuations that dominate the energy density at weak coupling correspond to the plasma 
particles and have momenta A; ~ T. For these "hard" fluctuations, 

/t ~ r, ((/?2)r ~ T^, 7r~5^- (54) 



Thus, at this scale, perturbation theory works as well as at zero temperature (with expansion 
parameter ~ y^, or rather a = /A-k). 
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The next "natural" scale, commonly referred to as the "soft scale" , corresponds to k ~ gT. 
We have then 

K ~ gT, {^'^)gT ~ gT'^, IgT ~ 9- (55) 

(In calculating {^'^)k for k <^ T, we have used the approximation ~ T/k, so that 
((/3^)K<cr ~ i^T.) Since ~ 5, perturbation theory can still be used to describe the 
self- interactions of the soft modes. However the perturbation theory is now an expansion 
in powers of g rather than g^: it is therefore less rapidly convergent. The emergence of this 
new expansion parameter is the origin of odd powers of g in the perturbative expansion of 
the pressure (such as the plasmon term ~ 5^). 

Another phenomenon occurs at the scale gT. While the expansion parameter 'jgT ~ g 
that controls the self-interactions of the soft fluctuations is small, the coupling between the 
soft modes and the thermal fluctuations at scale T is not: indeed g'^{ip^)T ~ (gT)'^, that 
is, the kinetic energy of the soft modes ~ (gT)'^ is comparable to their interaction energy 
resulting from their coupling to the hard modes, ~ g'^{(p^)T- Thus the dynamics of soft modes 
is non-perturbatively renormalized by their coupling to hard modes. This particular coupling 
is encompassed by the so-called "hard thermal loops" 66 1 . 



Finally, there is yet another scale, the "ultra-soft scale" k ~ g^T, at which perturbation 
theory completely breaks down. At this scale, we have indeed 

K ~ 5^r, {^'^)gT ~ g^T'^, Ig^T ~ 1- (56) 

Thus the ultra-soft fluctuations remain strongly coupled for arbitrarily small values of the 
coupling constant. This situation does not occur for a scalar field since a mass is generated 
at scale gT, which renders the contribution of the fluctuations at the scale g^T negligible. 
However this situation is met in QCD for the long wavelength, unscreened, magnetic fluctu- 
ations. 

These considerations suggest that the main difficulty with thermal perturbation theory is 
not so much related to the magnitude of the coupling constant (for the relevant temperatures 
it is not that large), but it is rather due to the interplay of degrees of freedom with various 
wavelengths, possibly involving collective modes. In some sense, field theories at high tem- 
perature are multiscale systems. At weak coupling the dynamically generated scales T, gT 
and g'^T are well separated. This allows for instance the organization of the calculation using 
effective field theory. However the scale separation disappears when the coupling is not too 
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small: then, the various degrees of freedom mix and the situation requires a different type of 
analysis. 

The exact renormalization group is ideally suited to cope with this type of situations. In 
particular, as we have seen, the BMW truncation scheme provides an excellent description 
of the momentum dependence of the 2-point function, from the low momenta of the criti- 
cal region, all the way up to the large momenta of the perturbative regime. One may then 
expect this method to capture accurately the contributions to the thermodynamical func- 
tions of thermal fluctuations from various momentum ranges, and hence handle properly the 
mixing between degrees of freedom that takes place as the coupling grows. Since it involves 
also non trivial momentum dependent self-energies, the method also encompasses implicitly 
effects related to the damping of quasiparticles, or their coupling to complex multi-particle 
configurations. 

Fig. [3] displays the pressure as the function of the coupling constant at the scale g{27rT). 
The various diverging curves labelled g^, - • • ,g^lng indicate the results of perturbative cal- 
culations, up to order g^ Ing (for recent high order calculations of the thermodynamics of the 



scalar field, see 



65|). These curves clearly illustrate the poor behavior of strict perturbation 



theory. The other curves correspond to various implementations of the exact renormaliza- 
tion group, as well as to a 2PI calculation based on a simple 2-loop skeleton [34]. The two 
curves labelled LPA correspond to two different choices of regulators: either the Litim regu- 
lator [6§], which is implemented only for three-momenta, with the (untruncated) sums over 
the Matsubara frequencies being performed analytically) [34], and an exponential regulator 
that affects both the momenta and the frequencies. The BMW approximation scheme is 
better justified when one uses an Euclidean symmetric four dimensional regulator, and the 
calculations reported here have been done with the exponential regulator (I42p . 

In contrast to the perturbative calculation, the calculations based on the renormalization 
group show a remarkable stability, and a smooth extrapolation towards strong coupling. As 
it turns out, the results obtained are not too different from those of the LPA, nor from the 
simple 2-loop 2PI approximation used in Ref. 3j]. In physical terms, both the LPA and 



the 2-loop 2PI approximation correspond to approximations where the degrees of freedom of 
the hot scalar plasma are massive quasiparticles. The new scheme goes beyond that simple 
picture. This stability of the results against improvements in the approximation suggests that 
the scheme that we are using to solve the NPRG equations may give already, at the level at 
which it is implemented here, an accurate representation of the exact pressure, and this over 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

g{2nT) 

FIG. 3: {Color online) Pressure as function of the coupling. The various resummation and RG 
methods compared are 2 PI, LPA with exponential regulator (exp), LPA with Litim regulator, and 
BMW. Shown are also perturbative results through order log g for the mass and g^ log g for the 
pressure. The g'^ and g^ logg curves for the pressure almost lie on top of each other, as do the BMW 
and LPA (exp) curves for the mass.. From Ref. [s^ 

a wide range of coupling constants. It also indicates that for such a system the quasiparticle 
picture is presumably robust. 



V. CONCLUSION 



The few applications that are presented in the previous section illustrate the power of 
the renormalization group, when coupled to an approximation scheme that allows for a 
determination of the full momentum dependence of the n-point functions. The equations 
that need to be solved are a priori complicated: these are flow equations which are at 
the same time partial differential equations (with partial derivatives with respect to the 
background field), with integral kernels that involve the solution itself. Still they can be 
solved at a rather modest numerical cost, using elementary numerical techniques. The 
studies presented here involved only the leading-order of a systematic approximation scheme. 
In the absence of any small parameter controlling the magnitude of successive orders, a 
study of the next order would be necessary in order to quantify the accuracy that has been 
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reached. However, the robustness of the leading order results can already be gauged from 
their weak residual dependence on the choice of the regulator. 
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